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There are many factors that contribute to the breaking of the spherical symmetry of a neutron 
star. Most notably is rotation, magnetic fields, and/or accretion of matter from companion stars. 
All these phenomena influence the macroscopic structures of neutron stars, but also impact their 
microscopic compositions. The purpose of this paper is to investigate the cooling of rotationally 
deformed, two-dimensional (2D) neutron stars in the framework of general relativity theory, with 
the ultimate goal of better understand the impact of 2D effects on the thermal evolution of such 
objects. The equations that govern the thermal evolution of rotating neutron stars are presented in 
this paper. The cooling of neutron stars with different frequencies is computed self-consistently by 
combining a fully general relativistic 2D rotation code with a general relativistic 2D cooling code. 
We show that rotation can significantly influence the thermal evolution of rotating neutron stars. 
Among the major new aspects are the appearances of hot spots on the poles, and an increase of the 
thermal coupling times between the core and the crust of rotating neutron stars. We show that this 
increase is independent of the microscopic properties of the stellar core, but depends only on the 
frequency of the star. 

PACS numbers: 04.40.Dg; 21.65.Cd; 26.60.-c; 97.60.Jd; 



I. INTRODUCTION 

The cooling of neutron stars has been used by many au- 
thors jTHHj as a way of probing the internal composition 
of these objects. Such studies rely on the fact that the 
physical quantities relevant for the cooling (specific heat, 
thermal conductivity, and neutrino emissions) strongly 
depend on the microscopic composition, so that differ- 
ent models lead to different thermal evolutions. For such 
studies, the predicted thermal evolution is compared to 
the observed data, with the ultimate goal of constrain- 
ing the microscopic properties of neutron stars [H [5J . An 
example of this is the recent analysis of the observed 
thermal behavior of the compact star in Cassiopeia A 
(Cas A) [TUl [HI [22] • In these studies the authors link 
the observed thermal behavior of this neutron star to the 
onset of superfluidity in the stellar core. Possible alterna- 
tive explanations have been suggested in [211 [25] , where 
the observed data was explained in terms of the nuclear 
medium cooling scenario [23] , and the a late onset of the 
Direct Urea process, triggered by the gravitational com- 
pression that accompanies a spinning-down neutron star 

[25]. 

In the standard approach, studies of the thermal evo- 
lution of neutron stars are performed by assuming that 
these stars are spherically symmetric. This assump- 
tion renders the thermal evolution calculations inherently 
one-dimensional (ID), which greatly simplifies the nu- 
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merical treatment. However, as pointed out by us in [25j . 
rotation may play a very important role for the cooling 
of neutron stars. In [23] the effects of rotation on the 
microscopic composition, and its consequences for cool- 
ing, were studied. In the work presented here we investi- 
gate the macroscopic aspects of rotation on the cooling of 
neutron stars in the framework of a fully self-consistent 
two-dimensional (2D) treatment. 

We shall consider the thermal evolution of rigidly ro- 
tating neutron stars. First we calculate the rotational 
structures of such stars, which is considerably more com- 
plicated than for spherically symmetric stars. The stellar 
structure is obtained by solving Einstein's field equation 
for a rotationally deformed 2D fluid (26] [27]. The nu- 
merical method used here is based on the KEH method 
[2HH30]- The equation of state used for computing the 
global neutron star properties and its composition is a 
relativistic non-linear mean field (RMF) model, whose 
parameters are adjusted to the properties of nuclear mat- 
ter at saturation density [51] . 

Once the structure of a rotating neutron star is com- 
puted, we solve the equations that govern the thermal 
evolution. The latter equations are re-derived for the 
metric of a 2D rotating neutron star. We show that the 
cooling of the object strongly depends on its frequency. 
At higher frequencies, rotating neutron stars show a sub- 
stantial temperature difference between the equator and 
the pole. We also show that the time scale for the ther- 
malization of such stars increases with its frequency. 

This paper is organized as follows. In section [n] we 
discuss the energy balance and transport equations for a 
rotating neutron star. In section [TTT] we show the results 
for the surface temperature evolution of rotating neutron 
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stars. The internal temperature evolution is discussed in 
section IV Finally, in section [V] we present our conclu- 
sions. 



II. 2D THERMAL EVOLUTION 

Given the complex nature of rotating neutron stars, 
where effects like frame dragging (Lense-Thirring ef- 
fect) and additional self-consistency conditions (e.g. mass 
shedding, which sets an absolute limit on rapid rotation) 
are in place, carrying out a 2D study of the thermal evo- 
lution of such objects is a challenging task. First steps 
toward this direction were made in J25J 1512H55] . We note 
that the 2D cooling of magnetized neutron stars has been 
previously studied in [35J [3B], where the 2D effects of 
the magnetic field in the crust were investigated, and 
first steps towards a consistent magneto-thermal evolu- 
tion simulation were taken. The cooling simulations per- 
formed in [3"51 [35] are based on a standard spherically 
symmetric metric, since this study focused primarily on 
the impact of magnetic field effects on cooling. Differ- 
ently from [35j, we are presenting full 2D calculations of 
the thermal evolution of rotating neutron stars, taking 
into account the 2D metric of a rotating fluid distribu- 
tion. The latter can be written as l2"51 
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where e 2 * = e 2i - a+ ^r 2 sin 2 9 and e 2 " = e 2 ^"^. The 
quantities v, <f> and uj denote metric functions, and N v ac- 
counts for the dragging of local inertial frames caused by 
the rotating fluid. All these quantities are functions of r 
and 9, and are implicitly dependent on N v . They need to 
be computed self-consistcntly from Einstein's field equa- 
tion, G 3 ' 5 = 87rT a ^, where T & P denotes the fluid's energy 
momentum tensor. 

The general relativistic equations of energy balance 
and transport are derived from the condition of energy- 
momentum conservation, and can be written as 
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where H = re^+^H^T, with H being the i-th com- 
ponent of the heat flux; T = e v T/T, with T being the 
temperature; k is the thermal conductivity; Cy is the 
specific heat; e is the neutrino emissivity; and the Lorentz 



factor T = (1 — C/ 2 ) -1 / 2 , where U is the the proper veloc- 
ity with respect to a zero angular momentum observer, 
given by U = (ft — N v )e^ . In the case of rigid body 
rotation, which is considered here, one has f2 =const so 
that Eqs. ^ to ^ reduce to 
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The standard cooling equations of spherically symmet- 
ric, non-rotating neutron stars are obtained from Eqs. 
^ through <HJ> for XI = and dgf = [26]. This work 
aims at solving Eqs. ^ to Q for the temperature distri- 
bution T(r, 9; t) of non-spherical, rotating neutron stars. 
The boundary conditions are obtained by defining H f at 
r = and R(9), and the heat flux Hg at 9 = 0,7r/2 and 
at r = R{9), with R(9) denoting the stellar radius. The 
star's initial temperature, T(r,9;t — 0), is typically cho- 
sen as T = 10 11 K. Equations ^ and Q can be solved 
for Hf and Hg, differentiated with respect to r and 9, 
respectively, and then substituted into Eq. ([2]), which 
leads to the following parabolic differential equation (for 
Vft = 0), 
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with the definitions r sin 6»e^+ 7 = and e- v+i = e""' 3 . 
Equation ([9| is solved numerically in combination with a 
fully general relativistic rotation code, which provides the 
necessary macroscopic input quantities. The microscopic 
input quantities (e.g. specific heat, thermal conductivity, 
neutrino emissivity) are obtained from the chosen model 
for equation of state. 

We should devote some time to discuss the microscopic 
model for the equation of state (EoS). In this work we 
used a non-linear relativistic mean-field model (param- 
eter set G300) [3T], which reproduces the properties of 
nuclear matter at saturation density. The underlying nu- 
clear lagrangian has the form [26j [27l [31] 

£ = ^2 V'sbvO^ - 9^ - gpP^l^f) - (m N - g a a)]ipB 
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FIG. 1: (Color online) Particle composition obtained from 
Eq. 1 10 1. The quantity p denotes the baryon number density, 



po is the number density at nuclear matter saturation density 
(0.16 fm~ 3 ). 



where B stands for protons (p), neutrons (n), and hy- 
perons (£, A, S). The interactions among these particles 
are described via the exchange of a, p, and u mesons. 
Their masses are m a — 550 MeV, m p — 769 MeV, and 
m u = 783 MeV and their coupling constants are given 
by g a = 9.1373, g p = 8.3029, and g u = 8.6324. The 
coupling constants of the self-interaction term of the a 
meson are b = 0.0033005 and c = 0.01529 [27]. The par- 
ticle population of neutron star matter computed from 
Eq. ( 10 ) is shown in Fig. [I] 

The G300 EoS leads to a relatively high proton frac- 
tion, triggering the direct Urea process (DU hereafter) 
|3"T] in neutron stars with masses above ~ 1.0 M Q . As 
discussed in previous studies [BTfTT] , the presence of the 
DU process enhances the cooling, leading to possible dis- 



agreement with observed data [4j [6] . This issue can be 
resolved by assuming that portions of the hadronic mat- 
ter in the cores of neutron stars are in a superfluid state 
[TUl EE] . The presence of an enhanced cooling process 
allows us to track the heat propagation more clearly in- 
side the rotating neutron star. This is due to the fact 
that, for rotating stars, the cold front originating from 
the stellar core is more intense (due to stronger neutrino 
emission, see [39]), and thus easier to track. 

As for the neutrino emission processes taking place in 
the core, we have considered the direct as well as the 
modified Urea processes together with bremsstrahlung 
processes. A detailed review of the emissivities of such 
processes can be found in reference [30]. In addition to 
the core, we also consider the standard processes that 
take place in the crust of a neutron star [40] . The spe- 
cific heat of the hadrons is given by the usual specific 
heat of fermions, as described in [3]. For the thermal 
conductivity we follow the calculations of |41) . 

In the present study, we are not considering the possi- 



bility that neutrons and protons may be in a superfluid 
and superconducting state, respectively [TO] [HJ [26] . 
The reason being that one obtains somewhat sharper 
temperature gradients inside the neutron stars if the neu- 
trino emission rates are not suppressed by superfluidity, 
which allows us to assess the effects of the 2D structure 
on the cooling of rotating neutron stars more clearly. The 
consequences of superfluidity /superconductivity, along 
with the Cooper-pair breaking/formation process, on the 
thermal evolution of rotating neutron stars will be pre- 
sented in a follow-up study. 



III. SURFACE TEMPERATURE EVOLUTION 

We now present the results for the thermal evolution of 
rotating neutron stars. We analyze the thermal evolution 
of three different stellar configurations, each one with the 
same central density but rotating at different selected 
frequencies, ranging from around 100 Hz to 800 Hz (see 
table 0. 

The choice of this frequency range is motivated by the- 
oretical calculations of the mapping between initial and 
final neutron star spins, likely spin-down mechanisms, 
and observational constraints from pulsars and supernova 
energetics. On this basis, it has been argued [32] that the 
initial pre-collapse central iron core periods inside of mas- 
sive stars, which give birth to neutron stars and pulsars, 
are on average greater than around 50 seconds. The as- 
sociated neutron stars would then be born with rotation 
periods greater than around 10 ms [42l [43] . Considering 
the many poorly understood features associated with the 
birth of rotating neutron stars, such as the estimated iron 
core spin rates, angular momentum profiles, progenitor 
masses, and general relativistic effects, the central iron 
core periods may be smaller by a factor of 5, however, in 
which case the birth periods of neutron stars increases to 
around 1000 Hz. Observationally, the most rapidly rotat- 
ing neutron star, PSR J1748-2446ad, rotates at 716 Hz 
[33], followed by the 642 Hz pulsar B1937+21 05]. Very 
recently, it has been noted that the low mass X-ray bi- 
nary SAX J1750. 8-2900 may have a spin frequency of 
601 Hz, and that this object may be the most luminous 
known neutron star |46| . 

Having said that, we now come back to the discussion 
of the properties of the neutron stars listed in table [T] The 
thermal cooling equation ^ was solved numerically for 
each star in table [I] and the results are shown in Figs. [2]- 
[4] where the redshifted temperatures at the poles and the 
equatorial belt of stars are plotted as a function of time. 

In Figure [2] we show the cooling of a slowly rotating 
(148 Hz) neutron star. At such low frequencies there is 
no significant difference between the polar and equatorial 
temperature evolution, indicating that such stars cool es- 
sentially like spherically symmetric stars (the eccentricity 
of this star is just 0.14, as shown in table [i]). The situa- 
tion is different for stars rotating at successively higher 



4 



TABLE I: Properties of neutron stars whose thermal evolu- 
tion is studied in this paper. All stars have a central den- 
sity of 350 MeV/fm 3 . Their properties are computed for the 
relativistic mean- field EoS G300 |31j . M denotes the gravi- 
tational mass, R e 



e = sj\ - (Rp/Re 
frequency, and Qk 



the equatorial radius 
the eccentricity. 



R p the polar radius, 



fi the star's rotational 
is the mass shedding frequency. 
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FIG. 2: (Color online) Redshifted temperature at star's pole 
and equator as a function of stellar age, for the 148 Hz star 
of table U 



frequencies, as shown in Figs. [3] and [4] One sees that as 
the frequency increases, so does the difference between 
the polar and equatorial temperatures, with the pole be- 
ing slightly warmer than the rest of the stellar surface 
due to the higher surface gravity there. 

It has been noted previously (35] [36] that anisotropic 
heat transport (caused by magnetic fields) also leads to 
the appearance of a hot spot in the magnetic pole. We 
note that in our case the hot spot is due to anisotropic 
heat transport due to rotation. Thus, if we were to in- 
clude magnetic field effects, we ought to obtain an even 
more pronounced hot spot at the pole (assuming that 
the magnetic and rotation axis are either the same, or 
differ only by a small inclination angle). A study which 
explores the impact of both effects combined is currently 
under way. 

Another noticeable difference concerns to the sharp 
temperature drop at ~ 100 years, which happens simul- 
taneously for the pole and equator of low frequency stars 
(i.e., spherical objects). This is not the case for stars ro- 
tating at higher frequencies, where a delay in the equato- 
rial temperature drop is observed, as shown in Figs. [3] and 
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FIG. 3: (Color online) Same as Fig. [2} but for the 488 Hz star 
of table |T] 
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FIG. 4: (Color online) Same as Fig. [5] but for the 755 Hz star 
of table □ 



[4] One can understand this sudden drop in surface tem- 
perature as the moment in time where the " cold front" , 
spreading from the core to the surface, has reached the 
stellar surface. (From this moment on the interior of the 
star is in thermal equilibrium.) If fast neutrino processes, 
like the direct Urea process, are active in the core, the 
cold front will be more distinctive, leading to a more pro- 
nounced drop in temperature. The time at which the cold 
front reaches the surface depends strongly on the prop- 
erties of the crust [39] . The cold front arrives first at the 
poles and then at the equator. This feature has its origin 
in the rotational deformation of the star: a flattening at 
the pole (leading to a thinner crust in this region) and 
an expansion at the equator (leading to a thicker crust). 
The thinner crust in the polar region, combined with the 
lower fluid velocity, leads to a more efficient heat prop- 
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FIG. 5: (Color online) Dependence of core-crust thermal cou- 
pling time, r, on rotational frequency of a neutron star. The 
quantity e c denotes the star's central density. Solid dots show 
the outcome of our numerical calculations, the lines show a 
4 th order polynomial fit to the data. 



agation, allowing the cold front to reach the pole more 
quickly, as shown in Figs. [3] and [4j This effect is less 
distinctive for stars with lower frequencies, since their 
shapes tend to be spherically symmetric and their crusts 
are more uniform. 

Keeping in mind that the core-crust thermal coupling 
time (t) depends primarily on the crust properties, and 
that the macroscopic properties of the crust are directly 
connected to the stellar rotation rate, we now investigate 
the dependence of r on the frequency of the star. In 
order to do that we calculate the cooling of two sets of 
stars, each set having central densities of 350 and 400 
MeV/fm 3 , and covering the range of frequencies from 
100 and 800 Hz. For each star we calculate t/tq, with 
To being the core-crust coupling time in the spherically 
symmetric case (i.e. Q = 0), and r the corresponding 
value for rotating neutron stars (defined as the moment 
that the cold front reaches the equator). The results are 
shown in Fig. [5} We find that t/to can be fitted well by a 
4 th order polynomial, and that the result is independent 
of the microscopic properties of the core, as can be seen 
m Fig. [5] The indicat ion is that this result should hold 
regardless of the EoS used for the core, as long as the 
crustal composition is the same. 

We note that in this study we have used the Baym- 
Pethick-Sutherland (BPS) model for the crust. The neu- 
trino emission processes taking place in the crust are 
Bremsstrahlung, e + e~ pair annihilation, and plasmon 
decay processes. A comprehensive overview of the neu- 
trino emission processes in NS can be found in [551 HO] • 
We have also assumed a non-magnetized neutron star 
atmosphere, composed of iron and iron-like elements. In 
which case the relationship between the mantle (T m ) and 
surface temperature (T s ) is given by the traditional rela- 



where T m g is in units of 10 8 K, T s q in 10 6 K, and g±4 is 
the surface gravity in units of 10 14 cm/s 2 . 



IV. INTERIOR TEMPERATURE EVOLUTION 

In order to better understand the thermal evolution of 
a rotating star, we now analyze the interior temperature 
evolution of such objects. The star rotating at ft = 755 
Hz will be used for this analysis, since its strong deforma- 
tion allows us to better evaluate the effects of rotation on 
cooling. The results of this analysis can be generalized 
to stars rotating at lower frequencies. 

Figure [6] shows a series of snapshots of the temper- 
ature profiles of the fl = 755 Hz star. Displayed are 
different cooling stages of this neutron star. For an early 
age of 0.1 year (Fig. [6] (a)) the core and crust are ther- 
mally decoupled, which can be clearly seen by the much 
higher temperature of the crust. Also clearly noticeably 
is the higher temperature of the pole, which is due to 
the greater surface gravity in this region. Fig. [6] (a) also 
shows that the core is much colder than the outer layers, 
which is a consequence of the fast core cooling via the 
DU process. The cold interior region will propagate like 
a "cold front" until it finally reaches the surface, at which 
point the star have become thermalized. The propaga- 
tion of the cold front can be observed in Fig.[6](b), which 
shows the temperature profile at a stellar age of 40 years. 



As already mentioned in section III the cold front prop- 
agates more efficiently in the polar direction, allowing 
this region to couple with the core more quickly than the 
equatorial regime (see also Figs. [2] to Eh. Figure [6] (c), 
(d) and (e) illustrate stellar epochs where the cold front 
is reaching the surface, first the polar region (Fig. [6] (c)), 
and then propagating downward towards the equatorial 
region (Fig.[6](d) and (e)). After 190 years, the stellar in- 
terior has reached thermal equilibrium, as shown in Fig. [6] 
(f). For ages greater than 190 years the pole is warmer 
than the equator, as was the case for Figs. [6] (a) and (b). 
The temperature difference is too small, however, to be 
noticeable in Fig. [6] (f). 

The surface temperatures shown in Fig. [6] are the red- 
shifted surface temperatures. This temperature is a func- 
tion of the mantle's temperature as given by Eq. (11) 

gang. 



V. DISCUSSION AND CONCLUSIONS 

The objective of this work was is to investigate in de- 
tail the cooling of rotationally deformed neutron stars 
in 2 dimensions (2D). Our investigation is fully self- 
consistent. Einstein's field equations of rotating compact 
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objects have been solved in combination with a 2D cool- 
ing code. Such studies arc extremely important if one 
wants to understand the thermal evolution of stellar sys- 
tems where spherical symmetry is broken. This might 
be the case for highly magnetized neutron stars, neutron 
stars that are spinning down from high to low rotation 
rates, as well as for accreting neutron stars which are 
being spun up to high rotation rates. This work rep- 
resents the first step towards a widely applicable, fully 
self-consistent treatment of the cooling of neutron stars 
whose spherical symmetries have been broken. As a first 
step in this direction, we considered here rigidly rotating, 
rotationally deformed neutron stars. The metric of such 
objects is significantly different from the metric of spheri- 
cally symmetric stars; we therefore re-derived the general 
relativistic equations that govern the thermal evolution 
of such stars. 

The microscopic model for the equation of state (EoS) 
which was used in our study is that of purely hadronic 
matter, with the whole baryon octet included. This 
model for the EoS allows for the direct Urea process in 
stars with masses above 1.0 M Q , which leads to enhanced 
stellar cooling. This choice for the equation of state was 
intentional, since the direct Urea process yields sharper 
temperature gradients inside the star, which allow us to 
track the energy transport more clearly. After obtaining 
a solid understanding of the thermal evolution of rotat- 
ing compact stars, our study will be applied in a future 
study to compact objects with more complex core compo- 
sitions, such as hadronic superfluids, boson condensates, 
or quark matter. 

Our results indicate that rotation (and hence the 2D 
structure) plays an important role for the thermal evolu- 



tion of neutron stars. The most important features can 
be summarized as follows: 

• The rcdshiftcd temperature of the pole is slightly 
higher than the temperature at the equator, which 
is due the difference in surface gravity between 
these two regions. 

• Thermal energy is more efficiently transported 
along the polar direction, due to the lower fluid 
velocities and thinner crust there. This allows a 
"cold front" , originating from the core, to arrive at 
the pole sooner than at the equator. 

With respect to the second item above, we have found 
that the increase in the core-crust coupling time is in- 
dependent of all stellar properties, except the frequency. 
This is consistent with the fact that the crusts of the 
neutron star at different rotational stages are very simi- 
lar regardless of their core composition, and that the only 
factor (in the present study) that alters the macroscopic 
properties of the crust is rotation. This is indicated by 
the increase in the core-crust coupling time as a function 
of frequency, which is very similar for objects with differ- 
ent central densities (and thus different particle compo- 
sitions). 
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FIG. 6: (Color online) Temperature profiles of a neutron star rotating at £1 — 755 Hz at different stages of stellar evolution, 
ranging from O.f (a), 40 (b), 60 (c), flO (d), f70 (e), to f90 (f) years. The colors represent the redshifted temperature. The 
stars are cut open to help visualize the interior temperature evolutions. 



